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ABSTRACT 

The Eastern Banded Structure (EBS) and Hydra I halo overdensity are very nearby (d ^ 10 kpc) 
objects discovered in SDSS data. Previous studies of the region have shown that EBS and Hydra I 
are spatially coincident, cold structures at the same distance, suggesting that Hydra I may be the 
EBS’s progenitor. We combine new wide-held DEGam imaging and MMT/Hectochelle spectroscopic 
observations of Hydra I with SDSS archival spectroscopic observations to quantify Hydra Ps present- 
day chemodynamical properties, and to infer whether it originated as a star cluster or dwarf galaxy. 

While previous work using shallow SDSS imaging assumed a standard old, metal-poor stellar popula¬ 
tion, our deeper DEGam imaging reveals that Hydra I has a thin, well-dehned main sequence turnoff 
of intermediate age (^5 — 6 Gyr) and metallicity ([Ee/H] = —0.9 dex). We measure statistically 
signihcant spreads in both the iron and alpha-element abundances of (J[Fe/H] = 0.13 ± 0.02 dex and 
^[a/Fe] = 0-09 ± 0.03 dex, respectively, and place upper limits on both the rotation and its proper 
motion. Hydra Ts intermediate age and [Ee/H] - as well as its low [g/Ec], apparent [Ee/H] spread, and 
present-day low luminosity - suggest that its progenitor was a dwarf galaxy, which has subsequently 
lost more than 99.99% of its stellar mass. 


1. INTRODUCTION 

The stellar tidal streams and substructure of the Milky 
Way (MW) provides an important observational window 
into the assembly of the Galaxy in a Universe dominated 
by dark energy and dark matter. In the standard AGDM 
model, numerical simulations show that massive galaxies 
like the Milky Way grow hierarchically, that is, from the 
continuous accretion and tidal destruction of smaller sys- 
tems such as dwarf galaxies and globular clusters (e.g.. 
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Despite the growing number of stream discoveries, only 
four stellar streams have known progenitors: the globu¬ 
lar clusters Pal 5, NGG 5466, NGG 288 and the Sagit- 
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tarius dwarf spheroidal. Developing a detailed picture 
of the origin of the Galaxy’s halo requires understand¬ 
ing the relative contribution of tidally-destroyed globu¬ 
lar clusters and dwarf galaxies. Searching for and quan¬ 
tifying the ages, metallicities, and kinematics of likely 
stream progenitors is therefore an important step to¬ 
wards a more complete picture of the hierarchical as¬ 
sembly of the MW’s halo in a cosmological context. 

The Eastern Banded Structure (EBS) stream (Eig- 
urefR is an ideal target for such a study, as it is relatively 
neaiTy {d ^ 10 kpc; Grillmair 2011) and therefor e may 
be close eno ugh for proper motion measurements. |Grill-| 
mair| ( |2006| ) discovered the EBS in SDSS imaging data 


as a relatively broad stell ar stream near t he Monoceros 
Ring stellar overdensity. Grillmair (2011) identified an 
extended double-lobed, ^ 2 sq. degree feature along 
the stream - called Hydra I - which has a suface den¬ 
sity larger than any other over density along the stream 
(white boxed region in EigurefB. Ba sed on its morphol¬ 
ogy and prominence, GrillmairpzOII) suggested that Hy¬ 
dra I may be the progenitor of the EBS stream. Using 
SDSS/SEGUE spectroscopic o bservations that are spa- 
tially coincident with Hydra I, Schlaufman et al. (2009) 
demonstrated the presence of a kinematically cold popu¬ 
lation of stars, with Uheiio ^ 85 km s“^ and the colors and 
magnitudes expected for an old, metal-poor population 
at the distance of the EBS. 

The origin and nature of the EBS/Hydra I features are 
complicated by their proximity to (and possible associ¬ 
ation with) a number of Gala ctic anti-center structures, 
including the Monocer os Ring ( Li et al.|2QI2 Slater et al. 
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Andromeda substructures (|Snettield et al. 2014 and ref¬ 
erences therein). The origin and relationshi p between 


these features remains intensely debated (see Xu et al 


2015 for a recent review). Eor example, it is currently un¬ 


clear whether the Monoceros Ring is tidal debris from an 
















































































Fig . 1.— Filtered and smoothed SDSS surface density map of the EBS stream and Galactic anti-center region (left', data from |Grillmair| 
[ 20 TTI 1 and the Hydra I overdensity (right). The Anti-Center Stream and Monoceros Ring are labeled for reference. The right panel is a 
zoom-in view of Hydra I (white box in the left panel) shown with the observing footprints of the DECam imaging (solid grey line) and 
MMT/Hectochelle spectroscopy (dashed grey line). 


accreted dwarf galaxy (Penarrubia et al. poo^ or has its 
origin in the Galactic disk - as either a normal warp/flare 
( Momany et al.|[2QQ6 ) or as the result of kicked-up stars 
from disk oscillati ons occurring from an encounter with a 
massive satellite dKazantzid is et al.||2QQ8 Price-Whelan 
et al. |2Q15| |Xu et al.| |2U15|). We consider our results in 


the context of these Gal actic anti-center structures in our 
discussion (Section [5^. 

Hydra I presents a unique opportunity to study a Milky 
Way satellite undergoing active disruption. In this work 
we determine Hydra Ps stellar population, chemical, and 
kinematic properties with the aim of exploring whether 
the object is a disrupting dwarf galaxy or star cluster. 
Previous studies of Hydra I region were limited by a 
small number candidate member stars (with large veloc¬ 
ity uncertainties), prohibiting detailed tests of Hydra Ps 
nature. In this work, we combine higher precision pho¬ 
tometric and spectroscopic data of Hydra I - from the 
Dark Energy Camera (DECam; Elaugher et al. 2Q15|) 


and M MT/Hectochelle spectrograph (Szentgyorgyi et at, 
2Q11[ ) - with with archival SDSS data to study Hydra I 


using a large sample of candidate stars. 

This paper is organized as follows. In Section I 
describe our observational data sets. In Section | 
discuss the selection of candidate member stars, and in 
Section we present the results of our analysis and the 
derived properties of Hydra I. We discuss possible scenar¬ 
ios concerning the nature of Hydra I in Section and 
present our conclusions in Section 


2. PHOTOMETRIC AND SPECTROSCOPIC DATA 
2.1. DECam Photometry 

Observations of Hydra I were obtained on 20-23 March 
2014 using the DECam imager installed on the CTIO 
Blanco 4m telescope. The DECam imager consists of 
sixty-two 2048x4096 imaging CCDs (chip gaps of ^ 50") 


arranged in a circular layout in the focal plane. On the 
Blanco telescope, this layout provides imaging over a 2.1 
degree diameter, resulting in a 3 square degree imaging 
area with 0.26" pixels. 

We imaged Hydra I using the ugri filters, but to fa¬ 
cilitate comparisons with other studies of Hydra I/EBS 
and the Galactic anti-center region, we only present the 
gr imaging in this work. We used a series of 6 x 180 sec 
exposures in each of the g and r filters and dithered be¬ 
tween exposures to fill in the chip gaps. Eigure shows 
the resulting footprint of the DECam imaging ( 9 deg^) 

on the Hydra I field. All imaging was obtained under 
photometric or near-photometric conditions. The mean 
EWHM of the image point-spread functions (PSEs) for 
the exposures ranged from 0.9" to 1.2". 

The data w ere reduced using t he NO AO Community 
Pipeline (CP; |Valdes et al.||2014|) and downloaded from 
the NO AO Science Archiv^ The CP performs both 
a standard instrumental calibration of the raw images 
(bias, overscan, cross-talk corrections; fringing correc¬ 
tions; dome flat and sky flat field corrections) and cor¬ 
rects for geometric distortions by reprojecting the images 
onto a common spatial scale. 

We performed photometry on each DECam exposure 
by separately analyzing each of the 62 reduced, repro¬ 
jected ctn 2 ^_^ing_^eDAOPHOTH/ ALLSTAR software 
suite ( Stetson|[l987[ 1994| ). We derive a catalog of point 
sources using cuts on the DAOPHOT CHI and SHARP 
parameters which vary as a function of magnitude depth. 
To ensure a robust characterization of sources, we require 
that object must be point-like in each of the gri images. 

The Hydra I/EBS region is located in the Sloan Digital 
Sky Survey (SDSS) footprint, and so we derive calibrated 
magnitudes by matching directly to SDSS Data Release 
8 (SDSS DR8;|Aihara et al.||20H|). We derived a photo¬ 
metric solution for each individual DECam exposure sep- 
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Fig. 2.— SDSS (left) and DECam (right) color-magnitude-number density (Hess) diagrams for stars within a 9 sq. deg region centered 
on Hydra 1. The diagrams have been slightly smoothed with a Gaussian filter to enhance low density features. The dashed white lines 
denote the CMP fi lter used to select candidates for spectroscopic follow-up with MMT/Hect ochelle. The iso chrone (solid white line, 
[Potter et al.||2008|> has an [Fe/H] = -1.8, age = 13 Gyr and shifted to a distance of 9.7 kpc ( |Grillm air||2011|). The median color and 
magnitude uncertainties are shown as white points. The higher precision DEGam photometry clearly shows tde detection of a blue MSTO 
at — r ~ 0.18 not apparent in the SDSS data. 


arately, after correcting the instrumental magnitudes for 
atmospheric extinction. We derived zero points and lin¬ 
ear color ter ms using the maximum l ikelihood approach 
described by Boettcher et ah (2013), using only point 
sources which have high b/N in both the DECam and 
SDSS imaging and SDSS colors of ^ — r < 1. We impose 
the latter constraint in order to improve the photometric 
calibration over the color region of interest for this study. 
Lastly, all calibrated apparent magnitudes presented in 


the reddening maps of Schlegel et al.j 

(19981) 

and the red- 

dening coefficients of bchlafly & F'ink 

beiner 

(2011 

). The 


Cam is E(B-V) 0.04 mag and varies by ±0.02 mag 
across the field of view. 

The final Hydra I DECam point source catalog was 
constructed by calculating the weighted mean magnitude 
of each point source from the individual measurements 
on each frame. We used the uncertainties in the cal¬ 
ibrated magnitudes (in flux space) as the weights and 
used an iterative sigma-clipping algorithm to reject any 
measurements which fell > 4a from the mean. The un¬ 
certainties in the weighted mean magnitudes were calcu¬ 
lated from the inverse square of the fiux errors. The final 
catalog contains ^ 17,000 sources between 17 < ^ < 22 
mag over a ^ 9 deg^ area. We present the entire photo¬ 
metric data set in Table which is available in machine 
readable format from the publisher. 

The number of point sources as a function of mag¬ 
nitude depth drops dramatically at ^ 23.5, r 22.8 

mag. We therefore estimate that our imaging is complete 
at ^ 1 mag brighter than these limits. A completeness 
limit at g ^ 22 — 22.5 is sufficient for the purposes of 
this study. As discussed in Section the main-sequence 
turn off (MSTO) in the color-magnitude diagram (CMD) 
of Hydra I is at ^ 18.5 (r 18.3) mag, significantly 

brighter than our estimated incompleteness depths. 


2.2. Hectochelle Spectroscopy 


We used MMT/Hectochelle to target stars throughout 
th e Hydra I overde nsity. Eollowing the previous analysis 
by Grillmair (2011|), we targeted stars from SDSS within 
0.1 rnag of a 13 Gyr, [Ee/H] = -1.8 i 


isochrone ( [Potter 


et al.|[2QQ8| ) at a di stance of 9.7 kpc (however, see our 


conclusions in § |43| regarding the age and metallicity of 
Hydra I). Eigurej^ shows the selection filter on the SDSS 
color-magnitude-number density (Hess) diagram of stars 
within the ^ 9 deg^ footprint of the DECam imaging. 
We chose a CMD filter with a relatively blue edge in or¬ 
der to explore the possibility of younger (bluer) stellar 
populations in Hydra I. Because the expected red giant 
branch (RGB) population of Hydra I overlaps with MW 
thick disk stars in color and magnitude (see Eigurej^, 
objects fainter than ^ = 18 mag were given priority for 
our Hectochelle observations in order to minimize con¬ 
tamination. 

Spectroscopic observations were obtained over 7 nights 
using the Hectochelle spectrograph on the 6.5m MMT. 
Hectochelle uses 240 fibers across a 1 degree diameter 
field of view to take single order spectra at a resolution 
of R ^ 38,000. We used the RV31 order selecting fil¬ 


ter to isolate the 150 A around the Mg b triplet (5167- 
5184 A). The field of view of Hectochelle is well suited to 
observing Hydra I. We obtained seven pointings around 
Hydra I, primarily in relatively bright or grey observ¬ 
ing time, with total exposure times ranging from 4800 
to 19,8000 sec per pointing. We targeted the areas on 
the primary overdensities and along the stream slightly 
southeast of the two main lobes. The footprint of the 
Hectochelle observations are shown in the right panel of 
Eigure We obtained spectra for a total of 796 stars 
with magnitudes in the range 17 < g < 21.3. 

The multifiber spectra were reduced in a uniform man¬ 
ner. Eor each field, the separate exposures were de- 
biased and flat fielded, and then compared before ex- 
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traction to allow identification and elimination of cos¬ 
mic rays through interpolation. Spectra were then ex¬ 
tracted, combined and wavelength calibrated using con¬ 
temporaneous exposures of a ThAr lamp (the spectra 
were not-linearized in dispersion). Each fiber has a dis¬ 
tinct wavelength dependence in throughput, which can 
be estimated using exposures of a continuum source or 
the twilight sky. A background sky estimate was made 
by combining data from fibers placed randomly in the 
focal plane, but excluding those which accidentally fell 
on a source. This sky spectrum was finely interpolated 
onto the object spectrum’s wavelength grid and then sub¬ 
tracted. 

We model th e Hectochelle s p ectra using the procedure 
introduced by Walker et al. (2015 hereafter W15) to 
model sky-subtracted Hectochelle spectra. The spectral 
model is given as: 


M{X)=Pi{X)t{X 


1 + 


Qm{^) T '^los 


( 1 ) 


where c is the speed of light and Pi{X) specifies the stellar 
continuum using an order-/ polynomial of the form 


Pi{X) =po+pi 


A-Ao 

A. 


+ P2 


r(A-Ao)l' 

A, 


+ .-.+PI 


r(A-Ao)l' 

A, 


( 2 ) 


The factor of T(A[1 + {QmW +'^ios)/c]) is a continuum- 
normalized template spectrum that is redshifted accord¬ 
ing to the line-of-sight velocity ('T’ios) and an order-m 
polynomial of the form 


Qm(^) = QO P Ql 


X-Xo 


A. 


q2 


(A - Ao) 


+ ... p Qm 


A. 

(A - Ao) 

A.S 


(3) 


that compensates for systematic differences between 
wavelength functions of target and template spectra (see 
W15 for details). Following W15, we choose m = 2, pro¬ 
viding sufficient flexibility to fit the observed continuum 
shape and to apply low-order corrections to the wave¬ 
length solution. We adopt scale parameters Aq = 5220 A 
and Ag = 60 A, such that —1 < (A — Ao)/A 5 < +1 over 
the entire range (5160 A - 5280 A ) considered in the fits. 

Again following W15, we generate template spectra T 
using the synthetic l ibrary that was use d for SSPP pa¬ 
rameter estimation (Lee et al. 2008a|b ). The SSPP li¬ 
brary contains rest-frame, continuum-normalized, stellar 
spectra computed over a grid of atmospheric parameters 
spanning 4000 < Teff/K < 10000 in effective temperature 
(with spacing ATeff/K = 250), 0 < \ogiQ[g/{cm/s‘^)] < 5 
in surface gravity (A logio[ 5 '/(cm/s^)] = 0.25 dex) and 
—5 < [Fe/H] < +1 in metallicity (A[Fe/H] =0.25 dex). 
The library spectra are calculated for a range in [a/Fe]; 
however, this ratio depends on metallicity. The library 
has [a/Fe] = -1-0.4 for spectra with [Fe/H] < —1, and 
then [a/Fe] decreases linearly as metallicity increases 
from — 1 < [Fe/H] < 0, with [a/Fe] = 0 for [Fe/H] > 0. 


The synthetic library spectra are calculated over the 
range 3000 — 10000 A at resolution 0.01 A /pixel, which 
we degrade to 0.05 A per pixel—twice as fine as our Hec¬ 
tochelle spectra. 

Our spectral model, M(A), is specified fully by a vec¬ 
tor of 13 free parameters. In order to fit the model and 
obtain estimates for each parameter, we follow W15’s 
Bayesian approach and adopt the same priors specified 
in their Table 2. We use the software pack age MultiNest 


(Feroz & Hobson 

to 

o 

o 

00 

Feroz et al.||2009) to sample the 

posterior probabi] 

[ity distribution functions (PDFs). As 


described in W15, the sampling of the PDFs allows one 
to evaluate the “Gaussianity” of the posterior PDFs by 
calculating the first four moments (mean, variance, skew, 
kurtosis) of the ID marginalized PDFs of each parame¬ 
ter. W15 note that in the high S/N regime, the poste¬ 
rior PDFs become increasingly Gaussian and therefore 
allow one to define a measure of quality control. That 
is, the Gaussianity measure provides confidence that the 
calculated variance on the velocity is a good measure of 
the Icr standard deviation (68% confidence interval). In 
practice, we keep only those spectra for which the model 
fits yield Gaussian shaped posterior PDFs as defined in 
W15 (see their Section 4.1 and Figure 4). 

While our fits to the Hectochelle spectra return simul¬ 
taneous estimates of stellar-atmospheric parameters Teff, 
log^, and [Fe/H], we And that these quantities exhibit 
systematic as well as random errors that are larger than 
those found, using the same technique on spectra from 
the same instrument, by W15. The reason for this dif¬ 
ference is that our current spectra were acquired in rel¬ 
atively bright conditions and our sky-subtracted spec¬ 
tra suffer significantly more residual contamination by 
scattered sunlight. Therefore, in this work we shall use 
only the velocities obtained from our Hectochelle spec¬ 
tra; for stellar-atmospheric parameters we will rely on 
the SDSS/SEGUE catalog. 

We examined the data quality and measurement reli¬ 
ability using both repeat observations and the pipeline 
data quality measures described above. To ensure a reli¬ 
able dataset, we adopted cuts on the radial velocity error 
(verr < 2.1 km s“^, or 3 times the median velocity error 
of 0.7 km s“^) and removed spectra with large residual 
contamination from scattered sunlight. After applying 
these cuts, 411 stars remained in the Hectochelle sam¬ 
ple. These data are presented in Table and include 
the DEGam magnitudes (corrected for Gmactic extinc¬ 
tion). The spatial positions of stars fainter than g = 18.5 
are shown with reject to the smoothed spatial map of 
Hydra I in Figure (left panel). The radial velocity er¬ 
ror of the Hectoch^e data as a function of magnitude is 
shown in Figure We obtain velocity errors less than 
^ 2 km s“^ to a limiting magnitude of ^ ^ 20.5, approx¬ 
imately two magnitudes below the main sequence turnoff 
(MSTO) of Hydra I. 


2.3. SDSS/SEGUE Spectroscopy 

We complemented our Hectochelle data with archival 
spectroscopy from SDSS, adopting the spectroscopic pa¬ 
rameters derived from the SEGUE Stell ar Parameter 
Pipeline (SSPP; Lee et al. 2008a, 2011). Unlike the 
MMT/Hectochelle target selection, the original SDSS 
survey spectroscopic selections were designed to meet a 
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Fig. 3.— Spatial positions of the Hectochelle/MMT sample (left panel) of targeted stars passing spectroscopic quality control cuts (see 
Section [ 2 ^ and the SDSS/SEGUE sample (right panel). Positions are shown on the smoothed surface density map of the Hydra I region 
(white boxed region in Figure Q. The white circle (r = 0.86 deg) shows the 2.3 sq. degree area centered on Hydra I corresponding to the 
DECam CMD shown in Figures andOnly stars fainter than g = 18.5 are shown for each sample. Stars are colored by their measured 
radial velocities. The white X shows the spatial center of Hydra I assumed in this study. 
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SDSS/SEGUE stars fainter than g = 18.5. 


SDSS/SEGUE 
• Hectochelle 



17 ^ 18 ^ 19 ^ 20 ^ 21 

g (mag) 


Fig. 4. — Radial velocity error vs g magnitude for the Hectochelle 
observations (black points) and archival SDSS/SEGUE data (grey 
points). Gross correlation failures resulting in large velocity errors 
are not shown. The dashed line indicates the velocity error cut 
(2.1 km s~ U ad opted in selecting the final, high-quality sample 
(see Section [ 2 ^ . 

number of scientific goals. We analyze the same SEGUE 
pointing as Schlaufman et al. (2009, 2011; B-8/PGI- 
9/PGII-21), but because the SDSS data are more het¬ 
erogeneously selected, we applied the same isochrone cut 
in color and magnitude used to select the Hectochelle 
targets - i.e., selecting only stars that might reason¬ 
ably be ^dra I members. These data are presented 
in Table 0 and include the DEGam magnitudes (cor¬ 
rected for^alactic extinction). Although SEGUE spec¬ 
troscopic stars are a biased sampling of stars within this 
color-magnitude selection, they should still yield a un¬ 
biased measurement of its kinematic properties and a 
lower limit on its spread in [Ee/H] and [a/Fe]. The right 
panel of Figure shows the spatial distribution of the 


3. PHOTOMETRIG AND SPEGTROSCOPIG SELEGTION 
OF HYDRA I STARS 

The DEGam photometry reveals a prominent MSTO 
feature in the Hess diagram (Figure at g 18.9, 
g — r 0.18 mag which is not apparent in the SDSS 
GMD. Because the DEGam imaging covers a relatively 
large area (compared to the spatial extent of Hydra I), 
we consider only a small region of 2.3 deg^ (radius of 
0.86 deg) centered on Hydra I in this study. This central 
region is shown on the smoothed spatial map of Hydra I 
in Figure Figure shows the DEGam GMDs of the 
central 2.3 sq. deg region and an equal-area region adja¬ 
cent to Hydra I for comparison. The comparison region 
was selected from the non-Hydra I DEGam imaging area 
located to the west of Hydra I. Only hints of the Hydra I 
MSTO and MS features which are present in the central 
pointing can be seen in the adjacent pointing, demon¬ 
strating that Hydra I is spatially concentrated. How¬ 
ever, because the adjacent comparison region is located 
relatively close to Hydra I, this region may not be truly 
sampling the background. We can therefore not rule out 
the possibility that the MSTO stars in the comparison 
region are simply EBS/Hydra I member stars located 
slightly outside the primary spatial overdensities. 

3.1. Minimizing Contamination in the Spectroscopic 
Samples 

The Hectochelle targeting employed a relatively wide, 
SDSS-based GMD filter, and given that we have clearly 
detected the MSTO of Hydra I, we can remove likely 
contaminants using their location on the higher preci¬ 
sion DEGam GMD. Stars with Hectochelle/MMT obser¬ 
vations are shown in Figure and are colored by their 
measured radial velocities. The SDSS/SEGUE sample is 
shown on the same GMD in Figure where stars are col¬ 
ored by their velocities ([^), iron ^undances ([^), and 
alpha abundances respectively. 
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Fig. 5. — DECam color-magnitude diagram for stars within the 2.3 sq. deg region (r < 0.86 deg) centered on Hydra I (top). The CMD 
of a equal-area region adjacent to Hydra I is shown for comparison (bottom). The area of the Hydra I CMD corresponds to the white circle 
shown in Fig[^ The Hydra I CMD shows a prominent MS with a relatively blue MSTO at ~ 18.9, — r ~ 0.18 mag. The Hydra 1 MS 

is less apparent in the adjacent field CMD, although the field does contain a significant number of stars similar to those at the MSTO in 
Hydra I. 


The radial velocities of stars along the MS are quite 
similar 90 km s“^; see Figs[6i and[^), whereas thick 
disk stars at brighter magnitmes show a larger veloc¬ 
ity spread and have higher iron abundances (—0.5 < 
[Fe/H] < 0). At fainter magnitudes {g > 18.5), the SDSS 
sample shows that both the iron abundance and velocity 
of MSTO stars are well correlated. In contrast, stars at 
slightly redder colors {g — r ^ 0.24) show a wider spread 


in velocity and may be on average more metal-poor than 
MSTO stars. 

We define a “faint sample” of Hectochelle stars using 
a narrower CMD filter shown by the blue dashed lines in 
Figure 1^. This filter has a red edge which is 0.06 mag 
bluer than our SDSS selection filter, minimizing contami¬ 
nation from redder stars {g — r^ 0.24). The filter also re¬ 
moves thick disk stars at brighter magnitudes {g < 18.5), 
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Fig. 6. — DECam CMD of stars within a 2.3 sq. deg region (r < 0.86 deg) centered on Hydra I (white circle in Fig[^. Panel (a) shows all 
stars (passing the quality control criteria) targeted by Hectochelle/MMT in central region, where the ^ints are colored by their measured 
velocity. The blue lines show the SDSS CMD filter used to select observing targ ets ( see also Figure 1^. Panel (b) shows the selection of 
the “faint sample” used for kinematic membership probabili ty an alysis (Sections |3.l) and |3.2| ). PaneT'(c) shows the derived membership 
probabilities Pm for the faint sample as described in Section [3^ Panel (d) shows the heliocentric velocity distribution of the bright (red 
histogram) and faint (black histogram) Hectochelle samples. 


slightly brighter than the MSTO. For our faint sample, 
we only include stars within the central 2.3 sq. degrees 
centered on Hydra I, resulting in a total of 122 stars. Fig¬ 
ure shows the velocity histogram of the Hectochelle 
faint sample. For comparison we show stars at brighter 
magnitudes {g < 18.5) which would have passed the nar¬ 
row CMD hlter. The contamination is readily appar¬ 
ent in the velocity histogram of bright stars, which only 
shows a small excess of stars at the measured systemic 
velocity of Hydra I {Vrad ^ 89 km s“^). By contrast, the 
faint sample shows that Hydra I is a signihcant overden- 
si ty of stars in veloci t y spa ce, consistent with the results 
of [Schlaufman et al.| (2009). 

I'he bDbb faint sample is selected without the use of 
a narrower, bluer CMD hlter due to the larger color 
measurement uncertainties. Although this may slightly 
increase contamination, we use the SSPP [Fe/H] mea¬ 
surements as an additional paramete r in o ur membership 
probability calculations (see Section 
stars fainter than g = 18.5 and witr 


3.2). We keep only 
im the 2.3 sq. de¬ 


gree region shown in Figure The SDSS faint sample 
contains a total of 42 stars. 

3.2. Spectroscopic Membership Probabilities 

We determine membership probabilities using the 
Expectation Maxim ization (EM) method described in 
Walker et al.|(|2009b|). The EM method uses ah available 


data to describe both the member and contamination dis¬ 
tributions. We assume a Gaussian velocity distribution, 
allowing the algorithm to return membership probabil¬ 
ities for each star along with the mean and variance of 
the member distribution. The spatial probability distri¬ 
bution is described as a monotonicahy decreasing func¬ 
tion of the stars’ distance from an assumed center be¬ 
tween the two primary lobes (RA=133.9°, Dec=3.6°; see 
Eig. |3| |Grihmair||20li|). The EM algorithm (and the de- 
riveomembership probabilities) implicitly assumes that 
there are no gradients in velocity, metahicity, velocity 
dispersion, or metahicity disp ersion. We use a Besan- 
con model (Robin et al. [2003) to describe the Galactic 
contamination in our EM algorithm. Eor the faint Hec- 
tochehe sample, we use only the spatial and velocity in¬ 
formation to determine membership probabilities. Eor 
the SDSS faint sample we also use the iron abundance as 
an additional parameter in determining memberships. 

The result of the EM algorithm applied to the Hec- 
tochehe faint and SDSS samples are shown in Eiguresj^ 
and ^ respectively. Each point is colored by its mem¬ 
bership probability (Pm)- The EM algorithm infers a 
systemic heliocentric velocity of V = 89.4 ±1.9 km 
s“^ and dispersion a = 8.4 ± 1.5 km s“^ for the Hec- 
tochehe faint sample, and Vsdss = 89.6 ± 5.3 km s“^ 
and CFSDSS = 9.1 ±3.1 km s“^ for the SDSS faint sam- 
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Fig. 7. — Identical to Figbut f or th e SDSS/SEGUE stars. The measured velocities, iron abundances, and alpha abundances from the 
SEGUE SSPP pipeline (see section |2.3| > are shown in panels (a), (b), and (c), respectively. The derived membership probabilities Pm (see 
Section are shown in panel (d). We discuss the mean values and spread in [Fe/H] and [a/Fe] for the high probability member stars in 
Section 
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The spatial positions of the Hectochelle and SDSS sam¬ 
ples, colored by membership probabilities, are shown on 
the smoothed surface density map of Hydra I in Fig¬ 
ure ^ The high probability members are well correlated 
with the photometric overdensities, although contamina¬ 
tion fro m likely non-mem bers is clearly present. As also 
noted by Grillmair (2011), the SDSS data show a clump 
of high Pm stars to the southeast of Hydra 1. Examining 
the larger EBS region in Eigureffl shows that these stars 
fall along the southeast extent of the EBS stream. 

4. RESULTS 

4.1. Limits on Rotation and Proper Motion 

The presence or lack of rotation in the Hydra I system 
is an important clue to its nature and possible ongoing 
destruction. To measure its rotation, we use a method 

( poo^lMo^ ) 

the Hectochelle 


simi lar to the one des c ribed by [Lane et al 
and 


Fig. 8. — Velocity versus radius for stars in the Hectochelle faint 
sample. Membership probabilities (Pm) were determined using 
velocity and spatial information with the EM method (Section |3.2| ). 
The mean velocity and velocity dispersion are listed. 

pie. Errors on the mean and dispersion wer e determined 
via bootstrapping. Schlaufman et al. (2009) found a sim¬ 
ilar mean velocity of 85 km s“^ from their analysis of the 
SEGUE pointing coincident with Hydra I and a higher 
velocity dispersion of 14.9 km s“^. 


Bellazzini et al. (2012) for stars m 
faint sample and SDSS saniples with membership proba¬ 
bilities Pm > 0.5. In this method, the sample is divided 
in two by a line that passes through the center of Hydra I, 
and a difference in the mean velocity of each subsample is 
taken. This velocity difference is calculated for dividing 
lines with position angles (PA) ranging from 0<PA<360 
deg. The analysis of both samples show no evidence for 
rotation at a statistically significant level. We find rota¬ 
tion amplitudes of ^ 3 —4 km s“^ with bootstrap-derived 
uncertainties of ^ 2 km s“^. We can therefore set a 3cr 
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Fig. 9. — Velocity versus radius (left panel) and [Fe/H] versus radius (right panel) for stars in the SPSS/SEG UE s ample. Membership 
probabilities (Pm) were determined using velocity, [Fe/H], and spatial information with the EM method (S ectio n |3.2| ). The mean velocity 
and velocity dispersions are listed. The mean and dispersion in the iron abundance are discussed in Section |4.2| 


Hectochelle stars 
g > 18.5 


4 

CD 


u 

CD 

Q 


• • X • 


136 


135 


134 133 

R.A. (deg) 


132 


131 


1.0 
0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
- 0.1 
lJo.o 





1.0 

0.9 

0.8 

0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0.0 


134 133 

R.A. (deg) 


131 


Fig. 10.— Spatial distribution of the spectroscopic samples (identical to Fign rej^ with stars colored by their EM membership probabilities 
(Pm)- Membership probabilities were determined as described in Section [3^ using velocity and spatial information; the SPSS sample uses 
[Fe/H] as an additional criteria for membership. The SPSS data show a clump of high Pm stars to the southeast of Hydra I overdensity. 
These stars fall along the extent of the EBS stream (see Figure]^ also |Grillmair|201 1| ) . 

.-1 


upper limit on the rotation of ^ 6 km s 


et al. 

20041 

tion o 

hlyd: 


Munn 


2008) to measure the ensemble proper mo- 
Hydra i stars in the Hectochelle sample. The typ¬ 
ical, per star, random measurement unce rtainty of this 
catalog is ^ 4 mas/yr (|Munn et ah 2004[). We find av¬ 


erage proper motions of cosd = —0.04 ± 0.57 mas/yr 
(2.4T34 km s“^) and JIs = —0.04T0.45 mas/yr (2.4T27 
km s“^) from Gaussian fits weighted by the EM mem¬ 
bership probabilities. Using the uncertainties in proper 
motions, we set a 3cr upper limit to the proper motion of 
/i ^ 2.2 mas/yr (132 km s“^). Although the measured 
proper motions are consistent with zero, they are very 
dif ferent than the p roper motion of Hydra I predicted 


by Grillmair (2011): /ic^cosd = —0.15 ± 0.57 mas/yr. 


Jls = —2.67T0.45 mas/yr for a prograde orbit. A combi¬ 
nation of increased numbers of stars and higher precision 
proper motions - perhaps from the Gaia mission - would 
be necessary for a refined proper motion measurement. 

4.2. [Fe/H] and [a/Fe] 

The mean iron- and alpha-abundances, and possible 
spreads, are important diagnostics for distinguishing be¬ 
tween a g lobular cluster or dwarf galaxy progenitor fo r 


the EBS (Willman fc Strader 2012 Kirby et al. 2013). 


As noted in Section|2.2| we use the SPSS sample to mea¬ 


sure the iron and alpha-element abundances of Hydra I 
stars, adopting the [Ee/H] and [a/Ee] values from SSPP 
and membership probabilities determined using the EM 
algorithm. 


Eigure 11 shows the [Ee/H] histograms for the SDSS 
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Fig. 11.— Iron abundance histograms of the SDSS data for the 
probable member stars (Pm > 0.5; black histogram) and probable 
non-member stars Pm < 0.5; red histogram). Using all stars, and 
weighting by membership probabilities, we find a mean [Fe/H] = 
—0.91 ± 0.03 dex and = 0.13 ± 0.02 dex when adopting the 

formal SSPP [Fe/H] uncertainties. Assum ing more conservative 
errors on [Fe/H] measurements (see Section 4.21 we find cr^pg/jj] = 
0.09 ± 0.03 dex. 


sample. Stars with Pm > 50% show a strong peak in 
iron abundance space, while lower probability stars show 
a wider spread in [Fe/H]. The alp ha-element \a/F e\ abun¬ 
dance measurements from SSPP (Lee et al. 2011|) provide 
an additional diagnostic on the star-by-star mernbership, 
and possible origin, of Hydra I. Figure 12 shows [tr/Fe] 
versus [Fe/H] for the SDSS sample. The^/Fe] distribu¬ 
tion of high probability member stars peaks near ^0.1 
dex, whereas the low membership probability stars show 
a broad, non-peaked distribution of [a/Fe] values. 

We determine the mean and dispersions of the iron 
and alpha-element abundances using the full sample of 
SDSS stars, adopting the EM membership probabilities 
as weights. We find mean abundances (with standard 
errors) of ([Fe/H]) = —0.91 ± 0.03 dex and ([a/Fe]) = 
0.14 ± 0.02. We measure statistically significant disper¬ 
sions in both abundances: cr[Fe/H] = 0.13 ± 0.02 dex and 
^[a/Fe] = 0.09T0.02 dex, respectively. We di scuss the im¬ 


plications of these measurements in Sections [5T] and [53 


How sensitive are the measured spreads in the abun¬ 
dances to the individual measurement uncertainties? To 
assess this, we recalculate the abundance spreads us¬ 
ing larger measurement uncertainties of 0.1 dex when 
the reported random errors on individual measurements 
are < 0.1 dex. Systematic comparisons to high reso¬ 
lution spectra and repeat SEGUE measurements have 
shown that the systematic errors on individual measure¬ 
ments are closer to ^ 0.1 dex, whereas the ran dom er¬ 
rors on high S/N s pectra are typically ^0.05 (Allende 


Prieto et al. 2008). With these conservative measure¬ 


ment uncertainties, we find cr[Fe/H] = 0.09 ±0.03 dex and 
^[a/Fe] = 0.08 ±0.02 dex. The use of larger uncertainties 
still yield statistically significant dispersions, and the val¬ 
ues themselves remain relatively unchanged at a spread 
of ^ 0.1 dex in the iron and alpha-elements. We note 
that adopting errors larger than ^0.12 dex effectively 
erases evidence of an abundance spread. 



-2.5 -2.0 -1.5 -1.0 -0.50 2 4 6 8 10 

[Fe/H] N 


Fig. 12.— [a/Fe] vs [Fe/H] for the SDSS sample. Black points 
indicate likely members of Hydra I (Pm > 0-5) and red points 
indicate likely non-members (Pm < 0.5). The error bars show the 
reported internal uncertainties in the measurements as reported in 
the SSPP catalog. Using all stars, and weighting by membership 
probabilities, we find a mean [a/Fe]^ 0.14 ± 0.02 and cr[Q,/[pg]] = 


0.09 ± 0.02 when adopting the formal SSPP [a/Fe] uncertainties. 
Assuming rnore conservative errors on the [a/Fe] measurements 
(see Section 4.2 I, we find cr[oi/Ye\ = 0.08 ± 0.02. 


4.3. Stellar Populations and Age Estimates 

Given the measurement of a mean iron abundance of 
[Ee/H] = -0.9 dex from the SDSS sample, we explore 
the age range of Hydra Fs stellar populations using the 
observed GMD of spectroscopic members and theoreti¬ 
cal isochrones. Although we used a color-rnagnitude cut 
to select spectroscopic targets (see Figures [^and|^, we 
are sensitive to MSTO stellar populations in Hydra I as 
young as 4 Gyr (for [Fe/H] = -0.9 dex). 

Figure shows the GMD of the Hectochelle and 
SDSS samples compared to isochrones with ages 5-6 Gyr, 
[Fe/H ] = -0.9 dex, and [a/Fe] = ±0.2 dex (jDotter et al. 
2008|). We infer a distance of d = 12.7±0.3kpc based on 
the prominent location of the MSTO in the GMD. We 
find that fitting the blue MSTO population of Hydra I 
stars {g — r ± 0.2 mag) requires a stellar population with 
an age of 5-6 Gyr, assuming a mean iron abundance of 
[Fe/H] = -0.9 dex. Ghanges of ±0.1 — 0.2 dex in the 
iron and alpha abundances result in older/younger ages 
by ^ 0.5 Gyr, which is below the age uncertainty of ^ 1 
Gyr implied by the photometric uncertainties in the DE- 
Gam g — r colors at the MSTO. The age of the ^5 — 6 
Gyr stellar populations correspond to a formation at red- 
shift ^ ^ 1 — 1.2 for the standard AGDM cosmological 


model ( 

Planck Gollaboration et al. 

2014). We discuss 

the implications of intermedia 
in Hydra I in Sections |5.1| and 
The bluest stars in Hydra I s! 

^e age stellar populations 

row good agreement with 


the ^5 — 6 Gyr isochrones and are clearly inconsistent 
with a single old stellar population (age > 7 Gyr). The 
Hectochelle sample shows a number of high probability 
stars at colors redder than the blue MSTO {g — r'> 0.2), 
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Fig. 13.— Comparison of the Hydra I CMD to theoretical isochrones. The top panel reproduces the CMD of the central region shown in 
Figure]^ The middle and bottom panels show the Hectochelle faint and SDSS samples, respectively, colored by the derived membership 
probab ffitie s Pm- For the SDSS sample, membership probabilities are derived using the [Fe/H] measurements as a second parameter (see 
Sectio n|3.2|) . The median uncertainty in the color is shown as a function of magnitude. Each panel shows fo ur isochrones for ages 5, 5.5, 6, 
and 7 Gyr for fixed values of [Fe/H] = -0.9, [Q;/Fe]=+0.2, at a distance of d = 12.7 kpc (|Dotter et al.|2008|). Hydra 1 contains a significant 
population of bluer stars consistent with ages of ~ 5 — 6 Gyr (see Section [T^ . 
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consistent with the expected colors of old, metal-poor 
MSTO stars. However, a number of factors limit our 
ability to draw conclusions about the presence of older 
stellar populations, including small numbers of candidate 
stars in both the Hectochelle and SDSS sample, as well 
as the limited/incomplete spectroscopic targeting of the 
reddest stars. We explored the possibility of a old metal- 
poor MSTO using a background subtracted Hess dia¬ 
gram of the region but found no statistically significant 
evidence for such a feature. Larger spectroscopic samples 
targeting the expect colors of older stars are necessary to 
further investigate the presence of complex stellar pop¬ 
ulations in Hydra I. We discuss the i mpli cations of the 
stellar populations analysis in Section [STT] 


4.4. Estimated Hydra I Luminosity 

We estimate a total of ^400 stars brighter than g = 
22 in Hydra I from a comparison of star counts in the 
overde nsity to the adjacent field (Figure 1^. As noted in 
Section [2Tj we expect no issues with com^eteness at this 
magnitude limit. We use the number count esti mate, in 
combination with a luminosity function (LF) from Hotter 
et al. (2008), to place limits on the total luminosity of 
Hydra i. I'he LF of Hydra I was created with a Salpeter 
IMF, an age of 5.5 Gyr, [Fe/H] = -0.9 dex, and [a/Fe] = 
0.2 dex. We calculated the total luminosity of Hydra I by 
normalizing the LF with our value of 400 stars brighter 
than g = 22. The normalized LF was then integrated to 
find the total luminosity of Hydra 1. We determine that 
the object has a total luminosity of ^ 1000 ± 200 Lq, or 
Mg —2.4±0.3 and Mr ^ —2.6±0.3. Our uncertainties 
are determined assuming Poisson statistics describe the 
number of star co unts. Using the SD SS-Johnson filter 
transformations of Jordi et al. (2006), we find My 
—2.5 ± 0.3 mag. 


4.5. Constraints on Dynamical Mass and Tidal Radius 

The measured velocity dispersion of Hydra I (a = 
8.4± 1.5), in combination with a half-light or other char¬ 
acteristic radius, can yi eld an estimate of the dynami¬ 
cal m ass of the system (Walker et al.||2009a Wolf et al. 
2010|). The spatial morphology of Hydra 1 indicates that 
the object is likely undergoing significant tidal distur¬ 
bance, so any measure of the dynamical mass is limited 
by the extent to which the present day stellar motions 
trace Hydra Fs underlying gravitational pote ntial. We 
discuss this assumption further in Section |5T| 

Within the inner r = 0.86 deg (190 pc) region of inter¬ 
est, we find a velocity dispersion a = 8.4 km s“^. If we 
assume this radius correspond to the half light radius of 
Hydra I, we can use Equation 2 from Wolf et al. (2010) 
to estimate the dynamical mass within the half light ra- 
dius, Mi/ 2 . We derive a dynamical mass estimate of 
^1/2 = 9 X 10^ Mq. If this mass estimate is within two 
orders of magnitude of the true dynamical mass, then 
Hydra I must have a significant dark matter com ponent 
given the luminosity of ^1000 Lq derived in |4.4| 

Given our stellar and dynamical mass estimates: Is 
the observation that Hydra I appears to be experiencing 
stellar mass loss consistent with the expected tidal radius 
of the object at its distance from the Galactic center? 
The tidal radius rudai of an object can b e approximated 
using Equation 2 from Bellazzini (2004): 


'^tidal 2 


^Hydra 


11/3 


2MMw(<^Hydra) 


d 


Hydra 5 


(4) 


where mnydra is the mass of Hydra I, dnydra is the 
galactocentric distance of Hydra I (^19 kpc), and 
Mmw {duydra) is the Milky Way mass within the dis¬ 
tance to Hydra 1. We calculate the Milky Way mass 
MMw{d}iydra) assuming an isothermal spher e with cir¬ 


cular velocity of Wire = 220 ± 40 km s ^ (Bellazzini 


20Mt 


We calculate the tidal radius (rudai) using both dy¬ 
namical and stellar mass estimates. Eor a dynamical 
mass of 9 X 10^ Mq, we find tidal radius of 1.7 deg 
(^ 370 pc). However, assuming a stellar mass to light 
ratio M^/L = 1 and L = 1000 Lq, we find a tidal radius 
of rtidai --0.1 deg (- 22 pc). 

The instantaneous tidal radii of both the “only stars” 
and the high dynamical mass scenarios appear to be con¬ 
sistent with our observations of Hydra 1. The small tidal 
radius, expected if Hydra I only contains stellar mass, 
is easily in good agreement with the spatial extent (^ 
1 deg radius) and asymmetric morphology of Hydra 1. 
The larger tidal radius of ^ 1 — 2 deg, expected from 
the estimated dynamical mass, is consistent with mass 
loss at Hydra Fs edges (possibly resulting in the EBS 
stream). It could also be consistent with significant tidal 
disturbance within Hydra Fs inner ^1 — 2 degrees, given 
that the instantaneous tidal radius may exceed its orbit 
averaged tidal radius and that we may be overestimating 
the dynamical mass. 

If, however, the system is unbound, we can estimate 
the dissolution time for the system based on the crossing 
time - the timescale for an unbound star to reach the 
observed extent of the object. For simplicity, we assume 
a characteristic size of 1 degree and that stars are mov¬ 
ing at their measured velocity dispersion (^ 8 km s“^). 
We find that the system would have a dissolution time 
of only ^ 30 X 10^ years. Given such a short dissolution 
timescale, it would be surprising to observe this stellar 
overdensity at the present day if Hydra I were truly un¬ 
bound. Although this provides general support for Hy¬ 
dra I as a bound object, additional kinematic data and 
numerical modeling would be necessary to draw more 
detailed conclusions. 


5. DISCUSSION 

We summarize the observed and derived properties 
from our study of Hydra I in Table In the following 
discussion, we consider the consistency of these proper¬ 
ties with different scenarios for the possible progenitor of 
the EBS stream. 


5.1. Is Hydra I a Disrupting Dwarf Galaxy? 

The hypothesis that Hydra I is a disrupting dwarf 
ga laxy has been con s idered previously, albeit indirectly, 

by ■ 

ha 


Schlaufman et al. ([2011|) in the discussion of their cold 
o objects (EGHUs) sample. They argued that the rel¬ 
atively high velocity dispersions of the EGHOs sample 
( (^vel) 20 km s ^) are difficult to explain with globu¬ 
lar cluster progenitors, w hich (on average) sh ow velocity 
dispersions of ^7 km s“^ (Ki mmig et al.|2Q15). However, 


see discussion below (Section |5.2p for a counterpoint to 
this broad dynamical argument. 
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TABLE 1 

Properties of Hydra I 


for an old, metal-poor stellar population (e.^. 


et al. 2013 
addition o; 


Brown et al. 2014 Weisz et al. 2014) 


Parameter 

Value 

Right Ascension 

133.9 deg 

Declination 

3.6 deg 

Galactic Longitude 

224.7 deg 

Galactic Latitude 

29.1 deg 

Distance 

12.7 ±0.3 kpc 

L 


- 1000 ± 200 Lq 

My 


~ —2.5 ± 0.3 mag 

Agf] 


~ 5 — 6 Gyr 

Hielio 1 


89.4 ± 1.4 km s-^ 

Vgsr^ 

g 

—55.5 ± 1.4 km s“^ 



8.4 ± 1.5 km s“^ 

Rotation AmplitudJ^ 

< 6 km s“^ 

([Fe/Hl>Pl 

-0.91 ± 0.03 dex 

f^[Fe/H 


0.13 ±0.02 dex 

{[a/Fe 


0.14 ±0.02 dex 

‘^[a/Fe 

P 

0.09 ± 0.03 dex 


Kirby 

this population to the observed metallicity 
distribution would significantly increase its dispersion in 
[Fe/H], making the value more consistent with other lu¬ 
minous dwarf galaxies. If Hydra I is indeed a dwarf 
galaxy, our non-detection of this population coul d b e 
partially due to our selection function (see Section 4.3). 
Another possibility is that Hydra I might have once pos¬ 
sessed such a population, but due to its larger spa tial 


extent (e.g., Battaglia et al. 2006 Bate et al. 2015), it 


^Y oung est stellar populations only (see Sec¬ 
tion 

^Derived from the faint Hectochelle sample. 

‘^Derived from the SDSS sample. 

If Hydra I is a dwarf galaxy, then the mass-metallicity 
and luminosity-metallicity relations for Local Group 
dwarf galaxies provide a means to estimate the stellar 
mass and luminosity of the stream progenitor (|Kirby 
et al.||2013). For a mean iron abundance of [Fe/H]= — 1 
dex, a dwarf galaxy stream progenitor would have been 
relatively massive - approximately ^ IO^Mq - sim¬ 
ilar to a present-day Fornax dSph. Given our stell ar es - 
timate of ^ 1000 Mq for Hydra I (see Sections |4.4|), 
this implies that a dwarf galaxy stream progenitor would 
have lost > 99.99% of its stellar mass. If Hydra I has un¬ 
dergone such substantial st ellar mass loss, it is unlikely 
that our analysis in Section [4^ yiel ds the true, gravita- 


has preferentially been tidally stripped compared to the 
more metal-rich subpopulation. We emphasize that our 
analysis does not rule out the possibility that Hydra I 
contains an old, metal-poor population. Additional spec¬ 
troscopic targeting of candidate stars in the expected re¬ 
gion of color-magnitude space and/or additional wide- 
held imaging would be necessary to hnd evidence for an 
older population. 

The presence of intermediate age st ellar population in 
Hydra I (ages ^5 — 6 Gyr; Section 4.3) is consistent 
with a dwarf galaxy progenitor hypothesis. Although 
MW dwarf spheroidals, on average, formed the majority 
of their stars prior to z ^ 2 (^ 10 Gyr ago), there is 
signihcant variation from galaxy to galaxy. Moreover, 
more massive dwarfs form a larger fract ion of their stellar 
mass at later times (jWeisz et al.||2014|). For example, ^ 
50% of the stellar mass in Lbrnax lorrned in the la st ^ 8 
Gyr and ^ 40% was formed in the last ^ 6 Gyr (Weisz 


tionally bound mass of the remnant. Smith et al. (2013) 


et al. 2014|). If Hydra I is a disrupting dwarf galaxy 
and was relatively massive at infall (as implied by the 
measured [Fe/H]), then the presence of an intermediate 
age stellar population is consistent with our expectations 
from the SFHs of Local Group dwarfs. 

5.2. Is Hydra I a Disrupting Star Cluster? 

The presence of ti dal streams associated w ith Milky 
Way GCs (e.g.. Pal 5;|Odenkirchen et al.|2001|) raises the 

n the hr 


|Grillmair| 

( 

Mill 

oled that, 1 

because 


discussed this possibil- 


find that the velocity dispersion of a dwarf galaxy rem- possibility that Hydra 1 could be a GC in the final throes 

nant is only a reasonable measure of the bound mass of destruction. . • 

when the object has retained > 10% of the initial mass. ity and conch 
When signihcant mass loss (> 98%) has taken place and 
only a small bound object remains - the remnant mass 
can be overestimated by as many as three orders of mag¬ 
nitude. We therefore don’t consider the large dynamical 
“mass” of Hydra I to provide support for a dwarf galaxy 
hypothesis. 

The iron abundance spread cr[Fe/H] is a potentially 
strong discriminator betw een dwarf galaxies and GGs 
(iWillman & Strader|2012|) . Gomparing a sample of Milky 
WT ^ 


Way dwarfs and GCs 


C^[Fe/H] 


illman & Strader (2012) hnd 
0.3 — 0.7 dex for dwarfs, while GCs less lu- 


on a highly eccentric orbit that brings it into the inner 
Galactic potential, dynamical heating from frequent en¬ 
counters with massive structures (such as giant molecu¬ 
lar clouds) or dark matter subhalos could be signihcant. 
Thus the large observed velocity dispersion for Hydra I 
could be due to dynamical heating rather than being ev¬ 
idence of a dwarf galaxy progenitor. 

To further explore the GG progenitor hypothesis, we 
consider the location of Hydra I in the age-metallicity 
relation (AMR) for MW globular and open clusters in 


minous than My = —10 have almost no measurable 
spread above the measurement uncertainties. Our es¬ 
timate of err Fe/H] — 0.1 ± 0.03 dcx is statistically sig¬ 
nihcant and consistent with a dwarf galaxy progeni¬ 
tor hypothesis. Although the brightest (i.e., classical; 
My < — 8) dwarfs generally have larger iron abundance 
spreads (^ 0.2 — 0.3 dex) compared to Hydra I, a rel¬ 
atively small value cr[Fe/H] is consistent with the appar¬ 
ent decre ase in er[Fe/H] with in creasing galaxy luminosity 
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( 2012 ). 


noted byfWillman & Strader 
All Milky Way dwarf galaxies with comprehensive 
studies of their stellar populations have shown evidence 


2011) and M31 ( Caldwell et al.||^v^......i;, w 

LIVIG and SMG (Fagel Taut vaisiene|| 1998). For the 
MW, GGs within 8 kpc of the Galactic center show a 
narrow range of old ages (^ 13 Gyr) but span a wide 
metallicity range. Outside a galactocentric distance of 
8 kpc, however, GGs exhibit an AMR more like that of 
the LMG and SMG, with a trend showing younger ages 
at higher metallicities. 

If Hydra I is a disrupting GG, then we expect it to con¬ 
tain (to first order) a single stellar population with the 
age of ^ 5 —6 Gyr as inferred from our stellar populations 
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Fig. 14.— Age vs [Fe/Hl for G 
anq/es; iDotter et al.|2010l |2011|) 

|Dias et al.||2Ui4|r Ages and [F ^ ^ 

clusters 4Y l'uc.~Whi tiug 1. and Pal 1 frouilHausen et al.T(l2013l>, 
|Valcheva et al.|(|2015|), andiRosenber g et al.lllUUtjalbl), res pectively, 
are also included and labeled. As in|uotter et al.f (|2011|), individ¬ 
ual GGs are identified by their galactocentric distance Rqc > 8 
kpc or Rgc ^ 8 kpc- The age range and mean metallicity of the 
stellar populations of Hydra I are shown as the grey boxed region. 
Very few known MW GGs or open clusters have ages and iron 
abundances similar to Hydra I. 


analysis (Section 4.3 
studied GCs in the in1 
These include Pal 1 ( 

). The MW has few known, well 
]ermediate age (^6 — 8 Cyr) range. 

Rosenberg et al. 1998a|b|), Pal 12, 

Terzan 7 dDotter et a 

l.||2010D, Whiting 1 (Valcheva et al. 

2015), and 47 Tuc (Hansen et al.|2013). With the excep- 


distances <18 kpc ( Harrisp996 ), similar to that of Hy¬ 


dra I (^Mw = 19 kpc). The youngest GCs of this subset 
(Pal 1 and Whiting 1), however, have higher iron abun¬ 
dances than Hydra 1. 

Outside of the MW, some evidence exists for massive 
star clusters with intermediate ages and metallicities. 
However, there are significant measurement uncertainties 
associated with these observations. For example, M31 
contains a large population of massive [Fe/H]^ —1 GCs, 
but age determinations from integrated light measure¬ 
ments (e.g.. Lick indices) at this metallicity are highly 
sensitive to the presence of blue horizontal branch stars, 
resulting in artificially young and thus high ly uncertain 
ages (see discussion inlCaldwell et al.||2011|). The Large 

TV /r n T /T TV /r y H k I I I I iiii . 


Magellanic Cloud (LMC) has a well "known “age-gap” in 


_ -ga 

its star cluster population (e.g., [Harris & Zaritsky|2009|) 
there are no massive clusters witn ages from ^ 5 — 12 Cyr 
and only one relatively low mass, intermediate metallic- 
ity cluster ( M ^ 10^ [Fe /H]^ —1) with an age 

^8 — 9 Cyr ( Mackey et al.|2QQ^ . In contrast, the Small 
Magellanic Cloud may contain a few M r\j 10^ -10^ Mo 
clusters c onsistent with th e age and metallicity of Hy¬ 
dra I (see Dias et aL]|2QlQ| for a recent review), although 
the wide range of age estimates for individual clusters 
highlights the large uncertainty in their observed ages. 

Given Hydra Fs outlying age and metallicity, relative 
to other known MW clusters, and the lack of strong ev¬ 
idence for analog clusters in M31 and the Magellanic 
Clouds, we conclude that the star cluster hypothesis is 


less likely than the dwarf galaxy hypothesis. If Hydra I 
is indeed a GC, the observed age and metallicity imply 
that it would have likely formed in dwarf galaxy (perhaps 
smaller than the LMC or SMC) which was subsequently 
accreted by the Milky Way. However, the lack of known 
GCs with ages and metallicities consistent with Hydra I 
places limits on the plausibility of such a scenario. 

5.3. Possible Association of Hydra I With a Larger 
Galactic Structure 

Lastly, we consider the possibility that Hydra I and 
EBS may be part of a larger Galactic structure - in par¬ 
ticular the Monoceros Ring stellar overdensity. The dis¬ 
tance and metallicity of Hydra I are consistent with pho¬ 
tometric and spec troscopic studies of the larger Mono - 
ceros Ring region (Conn et al.|2Q12 Meisner et al.|2Q12 ). 
Slater et al. (2014) recently presented a Pan-ST’ARRST 
stellar density map for the Monoceros Ring, constructed 
using MSTO stars from a ^ 9 Cyr population with 
[Fe/H] = -1 (colors in the range 0.2 < {g — r)o < 0.3). 
They argue that the EBS and Anticenter streams are ex¬ 
tended stellar arcs that are phy sically asso c iated with the 


Monoceros Ring. Conversely, Grillmair (2011) argued 


that the EBS stream is kinematically distinct from the 
Monoceros Ring, given the high eccentricity orbital solu¬ 
tion needed to explain EBS’s observed spatial curvature. 
The Monoceros Ring is believed to be on a nea rly circular 


orbit (Crane et al. 


2003 Martin et al. 2006). Although 


our study cannot shed additional light on whether either 
Hydra I or EBS might be associated with the Monoceros 
Ring, a more refined orbital solution for the EBS stream 
- using GAIA proper motions and a large sample of ra¬ 
dial velocity measurements - is a promising avenue for 
investigating this hypothesis. 

Given the proximity of Hydra I to the Monoceros Ring, 
it is possible we have isolated a younger population of 
stars recently speculated t o be part of this structure. 
Carballo-Bello et al. (2015) noted that along the line of 
sight to the IVlonoceros Ring (at I ^ 180°, 195°; b ^ 25°), 
there exists a population of stars bluer than the ^ 9 Cyr 
MSTO of Monoceros Ring stars. They suggested that 
this could be a ^ 4 Cyr population (with an iron abun¬ 
dance of [Ee /H] ^ —0.9 dex) asso ciated with the Mono¬ 


ceros Ring. |Carballo-Bello et al. (2015|) were unable to 
clearly detect a IVIST'O feature in the SDSS data but es¬ 
timated a distance of ^ 10 kpc for the population. The 
systemic velocities of the blue population - as well as 
Hydra I - are consiste nt with model predictions from 


Peharrubia et al. (2005) for the origin of the Monoceros 
Ring as an accreted dwarf galaxy. Although suggestive, 
additional precision photometry over a much larger field 
will be necessary to further explore any spat ial connec- 


(2015) detections 


tions between the [Carballo-Bello et al. 
and the EBS/Hydr 

The average [Ee/H] and [a/Ee] observed for Hydra I 


limits the possibility that it may just represent an over- 
densi ty of stars originating i n the MW’s halo or thick 


disk. Schlaufman et al. (2011) studied the SDSS/SEGUE 


observations of a set of cold halo substructures (ECHOs) 
one of which was Hydra 1. They compared the mean iron 
and alpha abundances of ECHOs to s amples of th i n disk , 
think disk, and halo stars compiled by Venn et al. (2004). 
They show, and we confirm, that their ECHO spatially 
associated with Hydra I is more iron-rich and less alpha- 
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enhanced than either the smooth halo or the thick disk 
in that direction ([Fe/H]haio ^ [<^/Fe]haio ^ 0.25]). 


6. CONCLUSIONS 

In this paper we present a wide-held DECam imag¬ 
ing and MMT/Hectochelle + archival SDSS spectro¬ 
scopic study of Hydra I - the double-lobed, ^ 2 deg^ 
halo overdensity previously suggested to be the progen¬ 
itor of the EBS stream - and quantify its present day 
chemo-dynamical properties. The results inferred us¬ 
ing both Hectochelle and SDSS/SEGUE data are con¬ 
sistent. Hydra I contains an intermediate age stellar 
population of ^ 5-6 Gyr with a mean iron abundance 
of [Ee/H] = —0.91 ± 0.02 dex and alpha abundance of 
[a/Fe] = 0.15 ± 0.03 dex. We hnd spreads in both the 
iron and alpha abundances of cr[Fe/H] — 0.10 ± 0.03 dex 
and cr[c^/Fej — 0.09 ± 0.03 dex, respectively. We conhrm 
previous observations that Hydra I is a kinematically cold 
structure, with a velocity dispersion of 8.4 ± 1.5 km s“^ 
at a systemic heliocentric velocity oiV = 89.4 ±1.4 km 
(Vgsr = —55.5 ± 1.4 km s“^). 

Hydra Fs stellar population appears initially inconsis¬ 
tent with earlier studies reporting that an old- and metal- 
poor stellar population yielded the strongest overd ensity 
in matched-hlter ing maps of the Hydra I region (Grill 
air||2006l |2011|). Old and metal-poor populations can 
have IVlbl'O populations as blue as those of intermediate 
age and [Ee/H] rich populations. Moreover, many ap¬ 
proaches to searching for and mapping streams in the 
MW’s halo assume that streams will be composed of 
old and metal-poor populations. Assuming a metal-poor 
population can thus lead to erroneous conclusions about 
a substructure’s age and distance. Complimentary spec¬ 
troscopy to measure the metal content of substructures 
provides a critical constraint on both their ages and dis¬ 
tances, and ultimately provides a more accurate picture 
of the accretion events that formed the Galactic halo. 

Hydra Fs intermediate age (^5 — 6 Gyr) and [Ee/H] 
(^-1 dex) make it an outlier among known MW globular 
cluster and open clusters. Conversely, Hydra Fs stellar 
populations, along with its spread in [Ee/H], is consis¬ 
tent with those of a Eornax-like dwarf galaxy progenitor. 
If Hydra Fs progenitor was a Eornax-like dwarf galaxy, 
then it must have lost > 99.99% of its stars to bear the 
low luminosity it displays today. Hydra I may thus be 
the first example of a system truly in its final throes 
of disruption, making it an excellent benchmark of tidal 
disruption models. 


Although evidence strongly points towards a dwarf 
galaxy interpretation of Hydra I, we cannot rule out the 
hypothesis that Hydra I originated as a star cluster or 
that it is physically associated with the Monoceros Ring. 
Higher S/N spectroscopic follow up of high probability 
member stars is a promising avenue for a more certain 
test of the dwarf galaxy hypothesis. With high resolu¬ 
tion spectra, the apparent [Ee/H] spread presented in this 
work could be measured with higher fidelity. Combined 
with Gaia proper motions and deeper wide-held imaging. 
Hydra I modeling will further test whether it could be 
kinematically associated with the Monoceros Ring and 
will constrain the impact of secular and environmental 
factors on its observable properties. 
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TABLE 2 

DECam Photometry 


R.A. 

(deg) 

Dec. 

(deg) 

go 

mag 

^90 

(mag) 

ro 

(mag) 

O'ro 

(mag) 

E(B - V) 
(mag) 

132.348949 

3.153806 

19.116 

0.015 

18.187 

0.009 

0.0318 

132.350640 

3.344006 

18.993 

0.013 

18.308 

0.009 

0.0339 

132.351339 

3.363441 

22.279 

0.042 

21.675 

0.033 

0.0340 

132.351990 

3.397592 

22.226 

0.068 

21.581 

0.038 

0.0340 

132.352598 

3.438227 

19.357 

0.012 

18.727 

0.009 

0.0338 

132.352717 

3.150096 

22.082 

0.037 

22.058 

0.050 

0.0319 

132.354884 

3.307660 

20.540 

0.014 

20.283 

0.013 

0.0341 

132.355516 

3.424209 

20.756 

0.016 

20.193 

0.011 

0.0339 

132.356115 

3.223843 

21.874 

0.027 

21.533 

0.025 

0.0334 

132.358492 

3.395449 

23.487 

0.106 

22.591 

0.070 

0.0341 


Note. — DECam photometric point-source catalog of the EBS/Hydra I 
region (see footprint in Figure^. The p hotometric data have been corrected 
for Galactic extinction (see Section [2m using the listed values oi E{B — V). 
We provide photometry only for stars with g — r < 1 due to the photomet¬ 
ric calibration method employed in this study. Point-source selection was 
performed using the DAOPHOT CHI and SHARP parameters. Table is 
published in its entirety in the electronic edition of The Astronomical Jour¬ 
nal. A portion is shown here for guidance regarding its form and content. 


TABLE 3 

Hectochelle/MMT Spectroscopy 


R.A. 

(deg) 

Dec. 

(deg) 

Vr 

(km s ^) 

(km s ^) 

90 

(mag) 

^90 

(mag) 

ro 

(mag) 

O'ro 

(mag) 

E{B - V) 
(mag) 

134.208389 

3.950312 

96.84 

0.44 

17.241 

0.005 

16.744 

0.003 

0.0499 

134.244110 

3.806066 

92.04 

0.53 

18.404 

0.003 

18.182 

0.002 

0.0479 

134.257034 

4.086879 

39.57 

1.10 

19.455 

0.005 

19.146 

0.004 

0.0512 

134.273636 

4.103878 

63.54 

0.44 

17.330 

0.005 

16.867 

0.003 

0.0510 

134.260468 

4.003685 

97.81 

0.44 

17.555 

0.005 

17.140 

0.003 

0.0505 

134.115387 

3.871296 

178.97 

0.35 

18.335 

0.004 

17.902 

0.003 

0.0514 

134.117737 

3.833826 

31.85 

1.47 

20.291 

0.006 

19.858 

0.004 

0.0501 

134.071915 

3.766812 

85.44 

0.97 

19.602 

0.004 

19.375 

0.003 

0.0486 

134.084824 

3.745721 

65.47 

0.34 

17.731 

0.003 

17.330 

0.003 

0.0475 

134.113266 

3.739799 

20.37 

0.46 

17.905 

0.004 

17.449 

0.003 

0.0464 


Note. — Hectochelle/MMT catalog of all 411 stars passing the spectroscopic quality control cuts 
and located within the SDSS CMD filte r sho wn in Figurej^. Only a subset of these spectra were used 
in our analysis of Hydra I (see Section [Sm. Helioc entr ic radial velocities {Vr) and Icr uncertainties 
were derived using the method des cribe d in Section |2.2| The DECam photometry has been corrected 
for Galactic extinction (see Section |2.l| ) using the listed values of E(B — V). Tablej^is published in its 
entirety in the electronic edition of The Astronomical Journal. A portion is shown here for guidance 
regarding its form and content. 


TABLE 4 

SDSS/SEGUE Spectroscopy 


R.A. 

(deg) 

Dec. 

(deg) 

Vr 

(km s-l) 

avu 

(km s-i) 

[Fe/H] 

(dex) 

0[Fe/H] 

(dex) 

[a/Fe] 

(dex) 

0[cx/Fe] 

(dex) 

90 

(mag) 

^90 

(mag) 

ro 

(mag) 

Oro 

(mag) 

E{B - V) 
(mag) 

132.572900 

3.247481 

54.42 

0.70 

-0.35 

0.051 

0.18 

0.014 

17.217 

0.006 

16.636 

0.006 

0.0348 

132.673340 

3.591954 

9.72 

0.80 

-0.20 

0.056 

0.13 

0.016 

16.880 

0.006 

16.377 

0.004 

0.0445 

132.678310 

2.828570 

269.38 

8.09 

-0.87 

0.066 

0.46 

0.018 

19.960 

0.013 

19.574 

0.009 

0.0304 

132.684040 

3.827853 

28.81 

1.49 

-0.47 

0.023 

0.31 

0.012 

18.131 

0.011 

17.584 

0.006 

0.0384 

132.685030 

3.294454 

44.80 

0.97 

-0.36 

0.040 

0.21 

0.014 

17.402 

0.006 

16.843 

0.003 

0.0343 

132.691000 

3.215672 

214.44 

2.48 

-1.57 

0.065 

0.24 

0.018 

18.027 

0.008 

17.526 

0.007 

0.0340 

132.703350 

2.779862 

26.07 

1.24 

-0.58 

0.031 

0.43 

0.013 

16.762 

0.010 

16.250 

0.007 

0.0307 

132.709140 

3.524834 

296.70 

1.97 

-1.01 

0.074 

0.34 

0.014 

18.124 

0.005 

17.600 

0.003 

0.0420 

132.734190 

3.252893 

53.05 

0.74 

-0.40 

0.011 

0.22 

0.013 

16.587 

0.005 

16.118 

0.004 

0.0339 


Note. — SDSS/SEGUE catalog of all 292 stars located in the Hydra I region and within the SDSS CMD filter shown in Figure |^. All 
spectroscopic parameters were taken from the SSPP pipeline. Only a subset of the se sp ectra were used in our analysis of Hydra I (see Section [3.1| ). 
The DECam photometry has been corrected for Galactic extinction (see Section [2m using the listed values of E{B — V). Table [^ is published 
in its entirety in the electronic edition of The Astronomical Journal. A portion is shown here for guidance regarding its form ana content. 














